Search Results

Multi-zone CFD simulations with detailed kinetics were used to model iso-octane HCCI experiments performed on a single-cylinder research engine. The modeling goals were to validate the method (multi-zone combustion modeling) and the reaction mechanism (LLNL 857 species iso-octane) by comparing model results to detailed exhaust speciation data, which was obtained with gas chromatography. The model is compared to experiments run at 1200 RPM and 1.35 bar boost pressure over an equivalence ratio range from 0.08 to 0.28. Fuel was introduced far upstream to ensure fuel and air homogeneity prior to entering the 13.8:1 compression ratio, shallow-bowl combustion chamber of this 4-stroke engine. The CFD grid incorporated a very detailed representation of the crevices, including the top-land ring crevice and head-gasket crevice. The ring crevice is resolved all the way into the ring pocket volume. The detailed grid was required to capture regions where emission species are formed and retained.

Meeting diesel engine emission standards for heavy-duty vehicles can be achieved by simultaneous injection of fuel and water. An injection system for instantaneous mixing of fuel and water in the combustion chamber has been developed by injecting water in a mixing passage located in the periphery of the fuel spray. The fuel spray is then entrained by water and hot air before it burns. The experimental work was carried out on a Rapid Compression Machine and on a Komatsu direct-injection heavy-duty diesel engine with a high pressure common rail fuel injection system. It was also supported by Computational Fluid Dynamics simulations of the injection and combustion processes in order to evaluate the effect of water vapor distribution on cylinder temperature and NOX formation. It has been concluded that when the water injection is appropriately timed, the combustion speed is slower and the cylinder temperature lower than in conventional diesel combustion.

Homogenous Charge Compression Ignition (HCCI) engines offer a good potential for achieving high fuel efficiency while virtually eliminating NOx and soot emissions from the exhaust. However, realizing the full fuel economy potential at the vehicle level depends on the size of the HCCI operating range. The usable HCCI range is determined by the knock limit on the upper end and the misfire limit at the lower end. Previously proven high sensitivity of the HCCI process to thermal conditions leads to a hypothesis that combustion chamber deposits (CCD) could directly affect HCCI combustion, and that insight about this effect can be helpful in expanding the low-load limit. A combustion chamber conditioning process was carried out in a single-cylinder gasoline-fueled engine with exhaust re-breathing to study CCD formation rates and their effect on combustion. Burn rates accelerated significantly over the forty hours of running under typical HCCI operating conditions.

Exhaust gas rebreathing is considered to be a practical enabler that could be used in HCCI production engines. Recent experimental work at the University of Michigan demonstrates that the combustion characteristics of an HCCI engine using large amounts of hot residual gas by rebreathing are very sensitive to engine thermal conditions. This computational study addresses HCCI engine operation with rebreathing, with emphasis on the effects of engine thermal conditions during transient periods. A 1-D cycle simulation with thermal networks is carried out under load and speed transitions. A knock integral auto-ignition model, a modified Woschni heat transfer model for HCCI engines and empirical correlations to define burn rate and combustion efficiency are incorporated into the engine cycle simulation model. The simulation results show very different engine behavior during the thermal transient periods compared with steady state.

Modern diesel engines operate under injection pressures varying from 30 to 200 MPa and employ combinations of very early and conventional injection timings to achieve partially homogeneous mixtures. The variety of injection and cylinder pressures results in droplet atomization under a wide range of Weber numbers. The high injection velocities lead to fast jet disintegration and secondary droplet atomization under shear and catastrophic breakup mechanisms. The primary atomization of the liquid jet is modeled considering the effects of both infinitesimal wave growth on the jet surface and jet turbulence. Modeling of the secondary atomization is based on a combination of a drop fragmentation analysis and a boundary layer stripping mechanism of the resulting fragments for high Weber numbers. The drop fragmentation process is predicted from instability considerations on the surface of the liquid drop.

One of the major problems limiting the accuracy of piezoelectric transducers for cylinder pressure measurements in an internal-combustion (IC) engine is the thermal shock. Thermal shock is generated from the temperature variation during the cycle. This temperature variation results in contraction and expansion of the diaphragm and consequently changes the force acting on the quartz in the pressure transducer. An empirical equation for compensation of the thermal shock error was derived from consideration of the diaphragm thermal deformation and actual pressure data. The deformation and the resulting pressure difference due to thermal shock are mainly a function of the change in surface temperature and the equation includes two model constants. In order to calibrate these two constants, the pressure inside the cylinder of a diesel engine was measured simultaneously using two types of pressure transducers, in addition to instantaneous wall temperature measurement.

Variable Compression Ratio (VCR) technology has long been recognized as a method of improving Spark Ignition (SI) engine fuel economy. The Pressure Reactive Piston (PRP) assembly features a two-piece piston, with a piston crown and separate piston skirt which enclose a spring set between them. The unique feature is that the upper piston reacts to the cylinder pressure, accommodating rapid engine load changes passively. This mechanism effectively limits the peak pressures at high loads without an additional control device, while allowing the engine to operate at high compression ratio during low load conditions. Dynamometer engine testing showed that Brake Specific Fuel Consumption (BSFC) improvement of the PRP over the conventional piston ranged from 8 to 18 % up to 70% load. Knock free full load operation was also achieved. The PRP equipped engine combustion is characterized by reverse motion of the piston crown near top dead center and higher thermal efficiency.

We have developed a methodology for analysis of Premixed Charge Compression Ignition (PCCI) engines that applies to conditions in which there is some stratification in the air-fuel distribution inside the cylinder at the time of combustion. The analysis methodology consists of two stages: first, a fluid mechanics code is used to determine temperature and equivalence ratio distributions as a function of crank angle, assuming motored conditions. The distribution information is then used for grouping the mass in the cylinder into a two-dimensional (temperature-equivalence ratio) array of zones. The zone information is then handed on to a detailed chemical kinetics model that calculates combustion, emissions and engine efficiency information. The methodology applies to situations where chemistry and fluid mechanics are weakly linked.

An experimental study has been carried out to provide qualitative and quantitative insight into gas to wall heat transfer in a gasoline fueled Homogeneous Charge Compression Ignition (HCCI) engine. Fast response thermocouples are embedded in the piston top and cylinder head surface to measure instantaneous wall temperature and heat flux. Heat flux measurements obtained at multiple locations show small spatial variations, thus confirming relative uniformity of in-cylinder conditions in a HCCI engine operating with premixed charge. Consequently, the spatially-averaged heat flux represents well the global heat transfer from the gas to the combustion chamber walls in the premixed HCCI engine, as confirmed through the gross heat release analysis. Heat flux measurements were used for assessing several existing heat transfer correlations. One of the most popular models, the Woschni expression, was shown to be inadequate for the HCCI engine.

Despite remarkable progress made over the past 30 years, automobiles continue to be a major source of hydrocarbon emissions. The objective of this study is to evaluate whether variable exhaust valve opening (EVO) and exhaust valve closing (EVC) can be used to reduce hydrocarbon emissions. An automotive gasoline engine was tested with different EVO and EVC timings under steady-state and start-up conditions. The first strategy that was evaluated uses early EVO with standard EVC. Although exhaust gas temperature is increased and catalyst light-off time is reduced, the rapid drop in cylinder temperature increases cylinder-out hydrocarbons to such a degree that a net increase in hydrocarbon emissions results. The second strategy that was evaluated uses early EVO with early EVC. Early EVO reduces catalyst light-off time by increasing exhaust gas temperature and early EVC keeps the hydrocarbon-rich exhaust gas from the piston crevice from leaving the cylinder.

The results of an experimental investigation of the flow in the near wake of a generic Sport Utility Vehicle (SUV) model are presented. The main goals of the study are to gain a better understanding of the external aerodynamics of SUVs, and to obtain a comprehensive experimental database that can be used as a benchmark to validate math-based CFD simulations for external aerodynamics. Data obtained in this study include the instantaneous and mean pressures, as well as mean velocities and turbulent quantities at various locations in the near wake. Mean pressure coefficients on the base of the SUV model vary from −0.23 to −0.1. The spectrum of the pressure coefficient fluctuation at the base of the model has a weak peak at a Strouhal number of 0.07. PIV measurements show a complex three-dimensional recirculation region behind the model of length approximately 1.2 times the width of the model.

This paper presents the utilization of alternative correlation functions in the Kriging method for generating surrogate models (metamodels) for the performance of the bearings in an internal combustion engine. Originally, in the Kriging method an anisotropic exponential covariance function is developed by selecting optimal correlation parameters through optimization. In this paper an alternative nonparametric isotropic covariance approach is employed instead for generating the correlation functions. In this manner the covariance for spatial data is evaluated in a more straightforward manner. The metamodels are developed based on results from a simulation solver computed at a limited number of sample points, which sample the design space.

This paper presents the development of surrogate models (metamodels) for evaluating the bearing performance in an internal combustion engine. The metamodels are employed for performing probabilistic analyses for the engine bearings. The metamodels are developed based on results from a simulation solver computed at a limited number of sample points, which sample the design space. An integrated system-level engine simulation model, consisting of a flexible crankshaft dynamics model and a flexible engine block model connected by a detailed hydrodynamic lubrication model, is employed in this paper for generating information necessary to construct the metamodels. An optimal symmetric latin hypercube algorithm is utilized for identifying the sampling points based on the number and the range of the variables that are considered to vary in the design space.

Adjusting the Residual Gas Fraction (RGF) by means of Variable Valve Actuation (VVA) is a strong candidate for controlling the ignition timing in Homogeneous Charge Compression Ignition (HCCI) engines. However, at high levels of residual gas fraction, insufficient mixing can lead to the presence of considerable temperature and composition variations. This paper extends previous modeling efforts to include the effect of RGF distribution on the onset of ignition and the rate of combustion using a multi-dimensional fluid mechanics code (KIVA-3V) sequentially with a multi-zone code with detailed chemical kinetics. KIVA-3V is used to simulate the gas exchange processes, while the multi-zone code computes the combustion event. It is shown that under certain conditions the effect of composition stratification is significant and cannot be captured by a single-zone model or a multi-zone model using only temperature zones.

This paper summarizes the first phase of an experimental effort focused on developing a comprehensive understanding of the in-cylinder air/fuel mixing and combustion processes in spark-ignition engines using laser-based fuel distribution and combustion measurements. As part of this first phase, a semi-quantitative, laser-induced fluorescence, fuel distribution measurement technique was developed and demonstrated. The calibration, correction, and image analysis processes associated with the technique were shown to be comparatively simple and effective (relative to other analytical and empirical methods). The error associated with the technique was shown to be 5 - 10 % under vapor phase conditions. This work was applied to a port fuel injected optical engine, which was designed for optical access through the piston and cylinder liner under firing conditions.

A method is presented here that detects aircraft tail surface icing that might normally be unobserved by the flight crew. Such icing can be detected through the action of highly computationally efficient signal processing of existing sensor signals using a so-called failure detection filter (FDF). The FDF creates a unique output signature permitting relatively early detection of tail surface icing. The FDF incorporates a stable state estimator from which the icing signature is created. This estimator is robust to analytical modeling errors or uncertainties, and to process noise (e.g. turbulence). Excellent performance of the method is demonstrated via simulation.

DKAD is an automated analysis tool for evaluating dynamic characteristics of accessory drives. Rotation response analysis predicts natural frequencies and effects of crankshaft excitation. Lateral response of each belt span shows the effect of pulley run-out and parametric excitation. DKAD systematically allows a user to define a design and its operating conditions and then performs a sequence of analysis to visualize the rotational and lateral responses. It also allows a user to quickly explore and assess alternative designs. Belt layout and associated parameters can be saved in templates for future reference.

A fast and accurate journal bearing elastohydrodynamic analysis is presented based on a finite difference formulation. The governing equations for the oil film pressure, stiffness and damping are solved using a finite difference approach. The oil film domain is discretized using a rectangular two-dimensional finite difference mesh. In this new formulation, it is not necessary to generate a global fluidity matrix similar to a finite element based solution. The finite difference equations are solved using a successive over relaxation (SOR) algorithm. The concept of “Influence Zone,” for computing the dynamic characteristics is introduced. The SOR algorithm and the “Influence Zone” concept significantly improve the computational efficiency without loss of accuracy. The new algorithms are validated with numerical results from the literature and their numerical efficiency is demonstrated.

A comprehensive formulation is presented for the dynamics of a rotating flexible crankshaft coupled with the dynamics of an engine block through a finite difference elastohydrodynamic main bearing lubrication algorithm. The coupling is based on detailed equilibrium conditions at the bearings. The component mode synthesis is employed for modeling the crankshaft and block dynamic behavior. A specialized algorithm for coupling the rigid and flexible body dynamics of the crankshaft within the framework of the component mode synthesis has been developed. A finite difference lubrication algorithm is used for computing the oil film elastohydrodynamic characteristics. A computationally accurate and efficient mapping algorithm has been developed for transferring information between a high - density computational grid for the elastohydrodynamic bearing solver and a low - density structural grid utilized in computing the crankshaft and block structural dynamic response.

Simultaneous measurements of the radial and the tangential components of velocity are obtained in a high-speed, direct-injection diesel engine typical of automotive applications. Results are presented for engine operation with fuel injection, but without combustion, for three different swirl ratios and four injection pressures. With the mean and fluctuating velocities, the r-θ plane shear stress and the mean flow gradients are obtained. Longitudinal and transverse length scales are also estimated via Taylor's hypothesis. The flow is shown to be sufficiently homogeneous and stationary to obtain meaningful length scale estimates. Concurrently, the flow and injection processes are simulated with KIVA-3V employing a RNG k-ε turbulence model. The measured turbulent kinetic energy k, r-θ plane mean strain rates ( 〈Srθ〉, 〈Srr〉, and 〈Sθθ〉 ), deviatoric turbulent stresses , and the r-θ plane turbulence production terms are compared directly to the simulated results.